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Abstract 

Capillary  pressure  versus  saturation  curves  for  drainage  of  a  wetting  phase  were  measured  for  several  gas  diffusion  layers  that  are  commonly 
used  in  polymer  electrolyte  membrane  fuel  cells.  The  technique  employed  can  measure  capillary  pressure  curves  for  both  the  total  pore  network 
and  the  pore  network  consisting  of  only  hydrophilic  pores.  This  enables  the  determination  of  capillary  pressure  curves  directly  relevant  to 
the  study  of  gas  diffusion  layer  flooding.  The  overall  distributions  compared  well  with  mercury  intrusion  data.  It  was  found  that  the  pore 
size  distribution  for  the  hydrophilic  pores  were  similar  in  shape  to  the  overall  distribution  for  standard  substrate  materials.  Materials  with  a 
microporous  layer  did  not  follow  this  trend  and  the  microporous  layer  was  found  to  be  completely  hydrophobic.  Due  to  their  similarity,  the 
overall  and  hydrophilic  capillary  pressure  curves  for  all  materials  could  be  correlated  using  a  single  Leverett  /-function.  The  results  were 
described  by  several  standard  capillary  models,  the  parameters  of  which  can  be  further  used  to  predict  the  relative  permeability  of  the  phases. 
©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Increasing  the  maximum  power  density  is  a  key  objec¬ 
tive  of  polymer  electrolyte  membrane  fuel  cell  (PEMFC) 
research.  Improvements  in  catalyst  activity,  accessibility  and 
utilization  have  led  to  a  dramatic  increase  in  reaction  rates  [  1  ] . 
The  corresponding  increase  in  water  generation,  however,  has 
tended  to  negate  these  improvements  due  to  flooding  of  the 
cathode  with  liquid  water  at  higher  current  densities.  The 
presence  of  excessive  liquid  water  in  the  cathode  can  debil¬ 
itate  the  cell  in  two  ways.  On  a  microscopic  scale,  liquid 
water  covers  catalyst  particles,  thereby  increasing  the  mass 
transfer  resistance  and  reducing  oxygen  access  to  the  catalyst 
sites.  Macroscopically,  liquid  water  may  fill  the  pore  network 
within  the  gas  diffusion  layer  (GDL)  support  and  decrease  the 
effective  diffusivity  of  oxygen  through  the  layer.  To  counter 


*  Corresponding  author.  Tel.:  +1  888  4567x3415;  fax:  +1  519  746  4979. 
E-mail  address:  mfowler@uwaterloo.ca  (M.W.  Fowler). 

0378-7753/$  -  see  front  matter  ©  2005  Elsevier  B.V.  All  rights  reserved, 
doi:  10. 1016/j.jpowsour.2005.05.086 


the  latter  effects  of  liquid  water,  GDLs  are  typically  coated 
with  a  non-wetting  polymer  such  as  polytetrafluoroethylene 
(PTFE)  to  create  hydrophobic  surfaces  and  pores  throughout 
the  GDL  that  are  thought  to  remain  free  of  water  and  facilitate 
transport  of  oxygen  to  the  catalyst  layer. 

Pore  scale  phenomena  associated  with  the  movement  of 
liquid  water  and  the  interplay  between  the  vapor  and  liquid 
phases  have  recently  received  considerable  attention.  Numer¬ 
ous  PEMFC  models  have  been  proposed  [2-20]  and  recent 
reviews  are  available  [21,22].  However,  analysis  of  these 
models  has  made  evident  the  scarcity  of  constitutive  rela¬ 
tions  that  describe  the  distribution  and  capillary  flow  of  liquid 
water  in  GDL  materials.  To  date,  very  limited  experimental 
information  has  been  reported  on  the  dependence  of  capillary 
pressure  (pc)  on  water  saturation  (,vw)  or  on  the  dependence 
of  water  relative  permeability  (krw)  on  saturation  in  different 
GDL  materials.  This  hinders  progress  in  model  validation. 

Several  studies  focusing  on  the  morphological  properties 
of  the  porous  GDLs  have  reported  Hg-Air  capillary  pres- 
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Nomenclature 

A 

area  (m2) 

BC 

Brooks-Corey  model,  given  by  Eq.  (10) 

D 

diameter  of  fibers  in  GDL  model 

J 

Leverett  J- function,  defined  by  Eq.  (16) 

k 

permeability  (m2) 

k. r 

relative  permeability 

m 

VG  model  fitting  parameter 

mA 

mass  of  dry  sample  (kg) 

ms 

mass  of  saturated  sample  (kg) 

n 

VG  model  fitting  parameter 

Pc 

capillary  pressure  (N  m-2) 

Pcb 

breakthrough  pressure  used  to  fit  BC  and  VG 
model  (N  m-2) 

RMS 

root  mean  squared  error: 

,  n  / 1  ^observed  predicted  2 

RMS  —  1  V  t  / (Pc  -/  c  ’ 

JXIVIO  n  /  J  A  /  observed 

1=1  v  Pc 

s 

saturation,  defined  by  Eq.  (1) 

s 

spacing  between  fibers  in  GDL  model 

VB 

bulk  volume  (m3) 

Vp 

pore  volume  (m3) 

VG 

van  Genuchten  model,  given  by  Eq.  (9) 

Greek  letters 

£ 

porosity 

9 

contact  angle  (rad) 

X 

BC  model  fitting  parameter 

P 

density  (kgm-3) 

G 

surface  tension  (N  m_  1 ) 

<p 

ratio  of  substrate  pore  volume  to  total 
pore  volume 

Superscripts 

BA 

SGL10BA 

BB 

SGL10BB 

eff 

effective 

S 

free  substrate  of  SGL10BB,  not  intruded 
by  MPL 

Subscripts 

B 

bulk 

G 

gas 

Hi 

hydrophilic 

nwp 

non-wetting  phase 

r 

residual 

T 

total 

wp 

wetting  phase 

sure  curves  or  pore  size  distributions,  but  none  have  been 
concerned  with  investigating  the  dependence  of  capillary 
pressure  on  saturation  for  the  water-air  system,  as  required 
for  modeling  water  transport.  Lee  et  al.  [23]  and  Kong  et  al. 
[24]  used  mercury  intrusion  porosimetry  (MIP)  to  examine 
the  consequences  of  various  microporous  layer  (MPL)  man¬ 


ufacturing  methods  on  the  morphology  of  the  MPL.  Jordan 
et  al.  [25]  and  Passalacqua  et  al.  [26]  experimented  with  dif¬ 
ferent  types  of  carbon  black  in  the  MPL.  Antolini  et  al.  [27] 
used  MIP  to  study  the  effect  of  PTFE  addition  on  the  total 
pore  volume  of  MPL.  Other  works  have  focused  specifically 
on  the  GDL  properties  [28-30],  but  were  aimed  at  qualita¬ 
tively  explaining  PEMFC  performance  on  the  basis  of  GDL 
morphological  features.  In  none  of  these  studies  has  a  sys¬ 
tematic  investigation  of  the  capillary  behavior  in  GDLs  been 
conducted.  All  of  the  research  has  been  carried  out  using 
mercury  intrusion  porosimetry  (MIP),  which  cannot  distin¬ 
guish  between  the  hydrophilic  and  hydrophobic  components 
of  GDL  porosity. 

In  the  absence  of  GDL  specific  pc-sw  data,  various 
assumptions  and  approximations  have  been  employed.  One 
common,  but  poorly  justified  approach  [2-8],  has  been  to  use 
a  polynomial  fit,  originally  obtained  by  Udell  [31],  of  capil¬ 
lary  pressure  data  measured  by  Leverett  for  water  imbibition 
in  water- wet  unconsolidated  sand  packs  [32].  Acknowledg¬ 
ing  the  lack  of  applicable  data  for  GDLs,  some  researchers 
have  arbitrarily  assumed  a  linear  dependence  between  capil¬ 
lary  pressure  and  saturation  [10,11,14],  A  few  attempts  have 
been  made  to  model  /\-sWp  relationships  specifically  for 
GDLs.  Divisek  et  al.  [17]  have  fitted  an  empirical  relation¬ 
ship  to  data  generated  from  a  model  developed  by  Ustohal 
et  al.  [33].  Pisani  et  al.  [18]  and  Natarajan  and  Nguyen 
[12,13]  adjusted  parameters  of  a  capillary  pressure  model  to 
match  observed  polarization  behavior.  Only  Weber  et  al.  [  16] 
expressly  attempted  to  empirically  fit  mercury  porosimetry 
data  for  use  in  their  models.  These  various  approaches  for  cal¬ 
culating  pc-sw  curves  have  produced  a  wide  range  of  results, 
leading  to  predictions  of  capillary  pressure  that  range  over 
several  orders  of  magnitude. 

The  heterogeneous  nature  of  GDL  materials  requires 
that  attention  be  given  to  the  coexistence  of  hydrophilic 
and  hydrophobic  pores  spaces  within  these  media.  Few 
researchers  have  addressed  this  point  [10,16,18]  and  only 
Weber  et  al.  [  1 6]  have  attempted  to  model  it.  The  latter  authors 
have  assumed  that  hydrophilic  and  hydrophobic  pore  spaces 
have  separate  size  distributions,  but  that  each  has  the  same 
shape  as  the  overall  distribution,  while  the  proportion  of  each 
pore  type  is  dictated  by  the  amount  of  PTFE  in  the  GDL  and 
its  expected  effect  on  hydrophobicity.  Although  this  approach 
is  more  realistic,  the  assumption  that  the  hydrophilic  and 
hydrophobic  pore  size  distributions  have  shapes  similar  to 
the  overall  pore  size  distribution  has  not  been  tested  against 
experimental  data.  Furthermore,  to  implement  this  approach, 
information  about  the  relative  proportions  of  hydrophilic  and 
hydrophobic  porosity  is  required.  In  this  respect,  only  slightly 
more  information  is  available  since  few  studies  have  inves¬ 
tigated  the  amount  of  hydrophilic  porosity  in  PTFE-treated 
GDLs  [28,34,35], 

The  objective  of  the  present  work  is  to  determine  the  pore 
size  distribution  and  void  fractions  associated  with  the  overall 
and  hydrophilic  porosity  of  various  commercially  available 
gas  diffusion  layers.  Such  a  breakdown  of  the  overall  porosity 
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distribution  of  a  GDL  into  its  hydrophobic  and  hydrophilic 
components  has  not  been  previously  reported.  Here,  the 
dependence  of  capillary  pressure  on  wetting  phase  satura¬ 
tion  along  a  desaturation  path  (drainage)  is  directly  measured 
using  the  method  of  “standard  porosimetry”  (MSP)  [36]  with 
octane-air  and  water-air  as  the  working  fluid  pairs.  Measure¬ 
ments  with  octane,  which  strongly  wets  both  hydrophobic  and 
hydrophilic  surfaces,  are  compared  to  MIP  measurements. 
Measurements  with  water  enable  the  determination  of  the 
pore  size  distribution  of  hydrophilic  porosity  of  interest  in 
fuel  cell  modeling.  The  effect  of  temperature  is  not  investi¬ 
gated  in  the  present  study  despite  the  documented  dependence 
of  capillary  pressure  curves  on  temperature  [37,38].  Conse¬ 
quently,  the  capillary  pressure  data  obtained  are  not  directly 
applicable  to  operating  fuel  cells.  Nonetheless,  these  data 
constitute  a  first  step  towards  understanding  the  capillary 
behavior  of  water-air-GDL  systems.  Fitting  of  the  capillary 
pressure  data  using  the  Brooks-Corey  and  van  Genuchten 
models  is  also  discussed.  Finally,  in  light  of  the  MSP  results, 
capillary  pressure  data  specific  to  the  hydrophobic  MPL  are 
obtained. 


2.  Materials  and  methods 

2.1.  GDL  materials 

A  representative  range  of  GDL  materials  was  selected  for 
investigation  in  this  work,  including  cloth,  felt  and  paper  with 
and  without  a  MPL.  The  materials  used  were  SGL  1 OBA  plain 
paper  and  SGL  10BB  paper  with  microporous  layer  (SGL 
Carbon  Group,  Short  Hills,  NJ),  Toray  TGPH-090  plain  paper 
(Toray  Corp.,  Tokyo,  Japan),  E-TEK  Cloth  ‘A’  (E-TEK  Inc., 
Somerset,  NJ),  as  well  as  Lyflex  484C  and  352C  felts  (Lydall, 
Manchester,  CT).  The  relevant  properties  of  each  material  are 
listed  in  Table  1 . 

2.2.  Mercury  intrusion  porosimetry 

Mercury  intrusion  porosimetry  (MIP)  is  widely  used  to 
measure  the  distribution  of  accessible  pore  volume  by  entry 
pressure  for  pores  between  1  nm  and  1  mm.  In  an  MIP  exper¬ 
iment,  the  non-wetting  phase  (mercury)  saturation  of  an  ini¬ 
tially  evacuated  sample  is  measured  as  the  capillary  pressure, 


Pc,Hs>-Air>  is  gradually  increased,  providing  the  capillary  pres¬ 
sure  curve  for  drainage  of  a  strongly  wetting  phase  [40].  Pore 
size  information  in  terms  of  an  equivalent  cylindrical  capil¬ 
lary  size  may  be  extracted  from  capillary  pressure  data  using 
the  Young-Laplace  equation: 

^  2(7 1  |g— Air  COS  0\  Ig — Air  ^  ^ 

Pc,Hg— Air 


where  r  is  the  pore  radius,  OHg-Air  the  surface  tension  of 
the  mercury-air  interface  and  0Hg-Air  the  contact  angle  of 
mercury  on  the  solid  surface  of  the  pore  wall.  In  general, 
inferring  the  pore  size  distribution  from  the  results  of  a 
MIP  experiment  is  fraught  with  difficulties  stemming  from 
pore  accessibility  limitations  [41].  Larger  pores  shielded  by 
smaller  ones  are  only  intruded  by  mercury  at  the  capillary 
pressure  corresponding  to  entry  of  mercury  in  the  smaller 
pores  and  their  volume  is  incorrectly  attributed  to  smaller 
pores.  In  chemically  heterogeneous  materials  such  as  a  GDL, 
the  assumption  that  mercury  exhibits  a  single,  well-defined 
contact  angle  on  all  solid  surfaces  is  also  challenged.  Par¬ 
tial  coverage  of  the  carbon  fibers  in  a  GDL  by  PTFE  creates 
non-uniformity  in  0Hg-Air,  meaning  that  two  pores  of  similar 
size,  but  with  different  wettability,  are  intruded  at  different 
pressures.  Uncertainty  regarding  the  value  of  @Hg-Air  thus 
confounds  the  conversion  of  capillary  pressure  measured  with 
one  fluid  (e.g.  Hg-Air)  to  equivalent  capillary  pressure  for 
another  fluid  pair  (e.g.,  water-air),  on  the  basis  of  the  follow¬ 
ing  equation: 


Pc, 2  = 


^2— Air  COS  02-Air 

- Pc  1 

CTi— Air  cos  9l—Air 


(2) 


in  which  oq-Air  and  (72-Air  are  the  surface  tensions  of  the 
two  liquids,  6>i_Air  and  02-Air  the  contact  of  angles  of  the  two 
liquids  on  the  material  and  pcq  the  experimentally  measured 
capillary  pressure. 

To  obtain  an  estimate  for  the  contact  angles  to  be  used  in 
Eq.  (2),  the  contact  angle  of  each  liquid  on  the  surface  of  each 
GDL  was  measured  by  image  analysis  using  a  video  contact 
angle  system  (AST  Products  2500XE,  Billerica,  MA).  The 
repeatability  of  these  measurements  was  ±1°.  The  observed 
contact  angles  were  then  corrected  for  the  effects  of  surface 
roughness  and  porosity  using  the  Cassie-Baxter  equation  for 


Table  1 

Physical  properties  of  GDL  materials  tested3 


Material 

Type 

e 

k  (m2  x  1012) 

pB  (gmL  ') 

Thickness  (|Jim) 

PTFE  (wt.%) 

SGL10BA 

Paper 

0.88 

18.0  [29] 

0.22 

380 

5 

SGL10BB 

Paper  w/MPL 

0.84 

0.33  [29] 

0.30 

420 

5  (in  paper) 

Toray  090 

Paper 

0.78 

8.3 

0.44 

190 

0 

E-Tek  Cloth  'A' 

Cloth 

- 

6.3b 

0.33 

350 

0 

Lyflex  484C 

Felt 

0.83 

- 

0.30 

417 

3.1 

Lyflex  352C 

Felt 

0.74 

- 

0.47 

268 

6.2 

a  Given  by  manufacturer  except  where  noted. 

b  Calculated  by  the  method  of  Happel  [39],  with  fiber  diameter  of  5  p.m  determined  by  fitting  model  to  Toray  090,  and  porosity  taken  as  0.75. 
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unit  cell 


Fig.  1.  Schematic  diagram  of  fiber  structure  used  to  determine  f\  and/2  in 
Eqs.  (5)  and  (6). 


porous  surfaces  [42]  modified  for  surface  roughness  [43]: 
cos  90bs  =  bfi  cos  0Eff  -  f2  (3) 


Table  2 

Observed  and  estimated  contact  angles  of  mercury  on  GDL  materials 


Material 

Observed  8  [°] 

Effective  8  [°] 

SGL10BA 

148 

108 

SGL10BB/MPLa 

148/153 

108/119 

Toray  090 

146 

112 

E-Tek  Cloth  'A' 

147 

114 

Lyfiex  484C 

147 

110 

Lyfiex  352C 

148 

114 

Graphite 

124  [44] 

- 

PTFE 

150  [43] 

- 

a  For  a  drop  deposited  on  MPL  side.  A  porosity  value  of  0.72  was  used  in 
Eq.  (3),  based  on  calculations  detailed  in  Section  3.2.6. 


f2 


s 

S+D 


(6) 


The  results  of  this  analysis  are  given  in  Table  2,  where  the  esti¬ 
mated  effective  contact  angles  are  also  compared  to  reported 
contact  angles  of  mercury  on  smooth  graphite  and  PTFE 
[43,44], 


where  f?obs  is  the  measured  contact  angle,  b  the  ratio  of  actual 
area  of  contact  between  the  drop  and  the  solid  portion  of  the 
surface  to  the  projected  area  (b  =  ;r/2  for  cylinders)  and/] 
and/2  the  fractions  of  the  GDL  surface  occupied  by  fiber  and 
void,  respectively.  The  value  of  #Eff  thus  obtained  is  a  rough 
estimate  of  the  contact  angle  of  the  liquid  on  the  chemically 
heterogeneous  fibers. 

To  estimate  the  values  of/i  and/2  for  the  GDLs,  the  fiber 
network  model  of  Nam  and  Kaviany  was  used  [10].  In  this 
model  the  GDL  is  composed  of  a  stack  of  interwoven  screens 
(Fig.  1).  Spacing  between  each  parallel  fiber  is  equal  in  both 
directions  and  this  is  also  the  spacing  between  the  layers.  For 
simplicity  the  fibers  are  allowed  to  intersect.  The  unit  cell 
in  this  model  has  dimensions  of  S  +  D  in  all  three  directions, 
giving  a  volume  of  (S  +  D)3 .  To  calculate  porosity,  the  fiber 
volume  is  subtracted  from  the  volume  of  the  unit  cell  and  the 
following  relationship  is  obtained: 

VunitCell  2Vpiber  +  ^Intersection 
S  =  - 

VunitCell 

(S+D)3-2fD2(S+D)+|D3 
(S  +  D)3 

(l)3  +  3(|)2  +  3(|)  +  l-f(|)-f  +  § 
(£)3+3(!)2  +  3(!)  +  l 

(4) 


Eq.  (4)  may  be  solved  for  S/D  given  e,  thus  enabling  the 
calculation  of  the  fractions  f\  and  />  from  the  following 
expressions: 


h 


D 

S+D 


1 


+  1 


(5) 


2.3.  Method  of  standard  porosimetry 

The  method  of  standard  porosimetry  (MSP)  offers  two 
distinct  advantages  over  MIP  [36].  Since  mercury  is  a  non¬ 
wetting  fluid  to  both  carbon  and  PTFE  surfaces,  an  MIP 
measurement  cannot  distinguish  between  hydrophilic  and 
hydrophobic  porosity.  MSP  can  be  performed  with  any  wet¬ 
ting  fluid  as  the  working  liquid.  Using  water  as  the  sat¬ 
urating  fluid  permits  one  to  determine  the  distribution  of 
hydrophilic  pores  only.  The  distribution  of  the  overall  poros¬ 
ity  can  be  obtained  by  using  a  strongly  wetting  fluid  such 
as  octane,  which  fully  wets  both  graphite  and  PTFE  sur¬ 
faces  uniformly.  The  second  advantage  of  this  technique  is 
that  the  /c-vwp  data  are  obtained  directly  for  the  fluid-solid 
system  of  interest.  Therefore,  the  results  are  not  subject 
to  the  uncertainties  in  contact  angle  introduced  by  the  use 
of  Eq.  (2).  To  date,  application  of  the  MSP  technique 
has  been  limited  to  the  characterization  of  cermet  materi¬ 
als  in  solid  oxide  fuel  cells  [45]  and  Nafion  membranes 
[46], 

The  method  of  standard  porosimetry  is  based  on  the  prin¬ 
ciple  of  capillary  equilibrium.  When  two  partially  saturated 
porous  bodies  are  in  contact,  the  system  moves  toward  an 
equilibrium  state  where  the  capillary  pressures  of  the  liq¬ 
uid  in  both  bodies  are  equal.  MSP  exploits  this  phenomenon 
by  placing  the  unknown  sample  in  capillary  contact  with  a 
standard  sample  having  a  known  pc-sw p  curve.  A  schematic 
diagram  of  the  steps  required  to  obtain  the  pc-sw p  curve  of 
an  unknown  sample  is  shown  in  Fig.  2.  When  the  system  is 
deemed  to  be  in  capillary  equilibrium,  the  mass  of  each  sam¬ 
ple  is  determined.  Knowledge  of  the  mass  of  each  sample 
permits  the  determination  of  its  wetting  phase  saturation  swp, 
defined  as  the  ratio  of  the  volume  of  liquid  in  the  sample  to 
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Fig.  2.  Procedure  for  obtaining  GDL  capillary  pressure  curves  using  MSP. 


the  total  void  volume  of  the  sample,  i.e., 


where  ms  is  the  mass  of  the  saturated  sample,  m,\  the  mass 
of  the  dry  sample,  pi  the  density  of  the  liquid,  s  the  porosity 
of  the  pore  network  being  measured  and  Vb  the  sample  bulk 
volume.  The  capillary  pressure,  pc,  in  the  standard  sample 
corresponding  to  its  saturation  can  then  be  determined  from 
its  known  capillary  pressure  curve  (Fig.  2a).  At  equilibrium 
this  is  also  equal  to  the  capillary  pressure  in  the  unknown 
sample  (Fig.  2b).  This  represents  a  direct  measurement  of 
the  capillary  pressure  in  the  unknown  sample  with  the  fluid 
of  interest.  This  process  can  then  be  repeated  over  a  range 
of  saturation  values  to  yield  the  capillary  pressure  curve  of 
the  unknown  sample  (Fig.  2c).  The  saturation  is  changed  by 
allowing  slow  evaporation  of  the  working  liquid  while  the 
standard  and  sample  are  in  capillary  contact.  The  capillary 
pressure  curve  thus  obtained  corresponds  to  drainage  of  the 
wetting  phase. 

One  drawback  of  MSP  is  its  inability  to  measure  residual 
saturations.  When  the  wetting  phase  only  moderately  wets 
the  porous  material  (such  as  water  in  the  hydrophilic  pores 
of  a  GDL)  there  is  the  possibility  that  rivulets  of  liquid  will 
detach  and  become  disconnected.  There  is  no  way  to  differ¬ 
entiate  the  mass  of  these  water  rivulets  from  the  mass  of  the 
connected  liquid,  since  the  evaporation  of  the  connected  and 
disconnected  liquid  occurs  simultaneously  during  the  desat¬ 
uration  phase  of  the  experiment.  As  a  result,  the  measured 
capillary  pressure  is  associated  with  the  total  saturation  in 
the  sample. 

2.4.  Experimental  procedure 

The  MSP  tests  were  run  on  a  Porotech  Automated  Stan¬ 
dard  Porosimeter  (Vaughan,  ON),  which  is  a  unique  apparatus 
for  the  automated  execution  of  the  MSP  procedure.  Each 
GDL  sample  tested  was  a  20  mm  diameter  circular  section 
cut  from  the  supplied  material.  The  standards  fitted  onto  the 
top  and  bottom  of  each  sample  were  fabricated  porous  disks 
of  proprietary  composition,  perfectly  wettable  by  both  octane 
and  water.  The  working  fluids  used  were  ACS  grade  octane 
(99.99%)  and  deionized  water  (>18.0  Mf2). 

One  of  the  most  important  steps  when  using  MSP  was  to 
ensure  that  the  samples  were  fully  saturated  at  the  start  of  the 
measurement.  To  accomplish  this,  the  sample  to  be  tested  was 
placed  in  a  sealed  chamber  that  was  subsequently  evacuated 


to  a  pressure  of  10  mTorr.  The  chamber  was  then  flooded  with 
the  working  liquid.  When  octane  was  the  working  fluid,  this 
procedure  was  sufficient  to  fully  saturate  the  sample,  since 
octane  spontaneously  imbibes  into  all  pores  due  to  its  high 
affinity  for  both  graphite  and  PTFE  surfaces.  When  water 
was  the  working  fluid,  an  extra  step  was  required  to  ensure 
that  all  hydrophilic  pores  had  been  filled.  After  removing  the 
sample  from  the  vacuum  chamber,  it  was  placed  in  water  at 
80  °C  for  several  hours.  The  elevated  temperature  reduced  the 
contact  angle  of  water  on  PTFE,  allowing  deeper  intrusion 
into  the  material  to  reach  any  inaccessible  hydrophilic  pores. 
The  sample  was  returned  to  room  temperature  before  testing. 

A  Quantachrome  Poremaster  (Boynton  Beach,  FL)  was 
used  for  MIP  testing.  Triply  distilled  ACS  grade  mercury 
(99.99%  purity)  was  used  to  obtain  the  pore  size  distribution 
of  the  samples.  Each  sample  with  mass  of  about  0.20  mg  was 
cut  into  about  20  small  20  mm  by  5  mm  rectangular  tabs  to 
fit  into  the  penetrometer  cell. 

3.  Results 

3.1.  Total  porosity 

The  overall  porosity  of  the  GDL  materials  (with  one 
exception)  was  provided  by  the  manufacturers,  whereas  the 
hydrophilic  porosity  was  not.  Table  3  compares  the  poros¬ 
ity  values  given  by  the  manufacturers  with  those  obtained  by 
MIP  and  MSP  in  this  study.  There  is  good  agreement  between 
the  manufacturers’  listed  porosities  and  those  determined 
experimentally  by  MIP  and  MSP  with  octane.  The  porosities 
measured  in  this  study  differ  from  the  manufacturers’  values 
by  not  more  than  0.04,  with  a  maximum  deviation  of  5.1%. 
The  fact  that  octane  reports  correct  values  for  total  porosity 
indicates  that  it  fully  wets  all  surfaces  and  the  assumption 


Table  3 

Total  and  hydrophilic  porosities 


Material 

Total  porosity,  sj 

Hydrophilic 
porosity,  £Hi 
MSP-water 

Manufacturer 

MIP 

MSP-octane 

SGL10BA 

0.88 

0.91 

0.86 

0.63 

SGL10BB 

0.84 

0.87 

0.81 

0.45 

Toray  090 

0.78 

0.79 

0.74 

0.63 

E-Tek  Cloth  'A' 

- 

0.72 

0.79 

0.74 

Lyflex  484C 

0.83 

0.86 

0.84 

0.81 

Lyflex  352C 

0.78 

0.74 

0.82 

0.36 
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of  zero  contact  angle  is  valid.  Conversely,  the  fact  that  the 
hydrophilic  porosity  is  less  than  the  total  porosity  in  materi¬ 
als  without  PTFE  (Toray  090  and  ETEK  Cloth  ‘A’)  indicates 
that  water  does  not  perfectly  wet  the  virgin  carbon  fibers. 

3.2.  Porosity  distributions 

In  this  section,  the  measured  capillary  pressure  data  are 
analyzed  in  terms  of  commonly  used  /;c-.vWp  models.  In 
the  following,  the  capillary  pressure  is  defined  in  terms  of 
the  difference  between  the  wetting  and  non-wetting  phase 
pressures  as: 

Pc  —  Pnwp  Pwp  (8) 


Inspection  of  the  MSP  results  for  the  two  Lyflex  samples 
and  SGL  10BB  indicates  the  presence  of  bimodal  pore  size 
distributions  in  these  materials.  These  results  are  understood 
by  considering  that  Lyflex  GDLs  are  heavily  laden  with  car¬ 
bon  powder,  which  introduces  structural  heterogeneity  in  the 
form  of  a  uniformly  distributed  microporosity,  whereas  the 
SGL  10BB  sample  has  a  microporous  layer,  which  is  both 
chemically  and  structurally  dissimilar  to  the  GDL  substrate. 

3.2.2.  Correlation  of  capillary  pressure  data 

One  of  the  main  goals  of  the  present  investigation  is  to 
provide  relevant  data  for  pc-swp  relationships  in  GDLs.  Sev¬ 
eral  models  are  commonly  used  to  fit  pc—sWp  curves.  One  is 
the  van  Genuchten  (VG)  model  [49]: 


3.2.1.  Capillary  pressure  curves:  total  pore  volume 

The  pc-sWp  curves  for  the  overall  pore  network  were 
obtained  using  MIP  and  MSP.  Octane  was  used  as  the  working 
fluid  for  the  MSP  measurements.  MIP  data  were  converted  to 
equivalent  octane-water  capillary  pressure  using  Eq.  (2)  and 
the  values  of  contact  angle  given  in  Table  2.  The  agreement 
between  the  capillary  pressure  curves  for  the  six  GDL  sam¬ 
ples  obtained  by  the  two  methods  is  quite  good,  as  evident 
from  Fig.  3.  Using  Eq.  (1),  the  mean  pore  radius  of  Toray 
090  was  calculated  from  the  MIP  data  as  9  p,m.  On  the  same 
material,  Park  et  al.  [47]  have  measured  a  mean  pore  radius  of 

1 1  p.m,  whereas  Mathias  et  al.  [48]  have  reported  a  value  of 

12  |xm  for  the  slightly  thinner  Toray  060  without  PFTE,  both 
results  obtained  by  capillary  flow  porometry.  The  agreement 
with  the  literature  values  for  Toray  090  and  with  MSP  results 
for  all  samples  (see  Fig.  3)  lends  support  to  the  effective  mer¬ 
cury  contact  angles  used  in  the  interpretation  of  MIP  data. 


-hst-  — 

where  pc b,  m  and  n  are  fitting  parameters.  In  physical  terms, 
pCb  corresponds  to  a  characteristic  capillary  pressure  (break¬ 
through  pressure)  associated  with  the  first  formation  of  a 
sample-spanning  cluster  of  pores  invaded  by  the  non-wetting 
phase  [41].  Another  common  model  is  the  Brooks-Corey 
(BC)  model  [50]: 


Pc  >  Pcb 


GO) 


In  both  models,  swp  is  identified  with  an  effective  saturation 
defined  as: 
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Fig.  3.  Comparison  of  capillary  pressure  curves  for  the  six  GDL  samples  obtained  using  MIP  and  MSP  with  octane  for  (a)  SGL  10BA  (b)  SGL  10BB  (c)  Toray 
090  (d)  E-Tek  Cloth  'A'  (e)  Lyflex  352  and  (f)  Lyflex  484  (Q:  MSP  (octane);  — :  MIP). 
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Table  4 

Model  fitting  parameters  for  total  pc-^wp  curves 


Material 

Brooks-Corey 

van  Genuchten 

Pcb  (bar) 

X 

RMS 

Pcb  (bar) 

m 

n 

RMS 

SGL10BA 

0.0562 

1.29 

0.25 

0.0802 

0.6489 

2.848 

0.34 

Toray  090 

0.1020 

1.59 

0.24 

0.1502 

0.7262 

3.652 

0.17 

E-Tek  Cloth  A 

0.0733 

1.08 

1.4 

0.0990 

0.5929 

2.457 

0.72 

where  ,vwpr  is  the  residual  saturation  and  .vwpt  the  total  satura¬ 
tion. 

When  applicable,  these  models  offer  the  possibility  of  esti¬ 
mating  the  dependence  of  wetting  and  non-wetting  phase 
relative  permeabilities  on  saturation.  For  example,  on  the 
basis  of  the  van  Genuchten  model,  the  relative  permeabili¬ 


ties  may  be  estimated  as  follows: 

‘™P  =  4'P2a  -  (i  -  si';")")1'2  <i2) 

*n> wp  =  (1  -  Swp),,3(l  -  (13) 

On  the  basis  of  the  Brooks-Corey  model,  the  equations  for 
relative  permeability  are: 

£rwp  =  4p2A  (14) 

£rnwp  =  (l-Swp)2(l-4p2A)  (15) 

A  drawback  of  both  these  models  is  that  they  can  only  be 


applied  to  unimodal  pc-Swp  curves,  and  thus  are  not  appro¬ 
priate  for  the  SGL  10BB  and  the  Lyflex  felt  GDLs  (Fig.  3). 
Carbon  powder  is  incorporated  in  these  GDLs,  leading  to  a 


bimodal  pore  size  distribution.  Consequently,  capillary  pres¬ 
sure  curves  were  fit  according  to  the  BC  and  VG  models 
only  for  the  remaining  three  samples.  The  fitting  parame¬ 
ters  so  obtained  are  given  in  Table  4.  In  the  case  of  the  BC 
model,  data  points  below  pc=Pcb  (sWp  greater  than  ca.  0.8) 
were  excluded  since  this  model  is  not  intended  to  describe 
the  behavior  in  these  regions. 

A  comparison  of  the  fitted  curves  and  the  experimental 
data  for  these  samples  is  presented  in  Fig.  4.  As  expected, 
the  BC  model  does  not  agree  well  with  the  experimental 
data  at  high  wetting  phase  saturations.  The  VG  model  is  the 
better  choice  overall,  as  it  captures  the  smooth  rise  in  cap¬ 
illary  pressure  at  high  wetting  phase  saturations.  This  rise 
is  due  to  the  fact  that  a  large  proportion  of  the  pore  vol¬ 
ume  is  accessible  directly  from  the  external  GDL  surfaces, 
allowing  these  pores  to  be  penetrated  by  the  fluid  at  pres¬ 
sures  below  the  breakthrough  pressure  of  the  medium.  In 
hydrogeology  and  reservoir  engineering  applications,  where 
the  capillary  behavior  of  a  very  large  volume  of  materials  is 
of  interest,  this  part  of  an  experimental  capillary  pressure 
curve  is  ignored  since  it  is  not  representative  of  an  infi¬ 
nite  medium.  However,  this  region  is  an  intrinsic  aspect  of 


Water-Air  Capillary  Pressure  [Bar] 


Fig.  4.  Comparison  of  model  fit  and  experimental  capillary  pressure  curves  determined  by  MSP  with  octane  for  (a)  E-TEK  Cloth  ‘A’  (b)  Toray  090  (c)  SGL 
10BA.  Solid  lines  are  VG  model  and  dashed  lines  are  BC  model. 
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Fig.  5.  Comparison  of  octane-air  and  water-air  capillary  pressure  curves  obtained  using  the  MSP  technique  for  (a)  SGL  10BA  (b)  SGL  10BB  (c)  Toray  090 
(d)  E-Tek  Cloth  ‘A’  (e)  Lyflex  352  and  (f)  Lyflex  484  (•:  hydrophilic  pore  network  (water);  Q:  total  pore  network  (octane)). 


GDL  behavior,  particularly  at  the  GDL-gas  channel  inter¬ 
face,  and  so  should  not  be  ignored.  Considerable  error  asso¬ 
ciated  with  the  fitting  of  the  VG  model  is  mostly  observed 
in  the  high  drwp  region.  If  only  wetting  phase  saturations 
below  95%  are  considered,  the  error  is  less  than  10%  for  all 
samples. 

3.2.3.  Capillary  pressure  cun’es:  hydrophilic  pore 
volume 

The  capillary  pressure  curves  associated  with  the  drainage 
of  water  from  the  network  of  hydrophilic  pores  were  obtained 
using  MSP  with  water  as  the  working  fluid.  Fig.  5  shows  the 
pc—s nwp  curves  for  the  hydrophilic  pore  network  for  each 
GDL  under  study.  The  capillary  pressure  curves  for  drainage 
of  octane  from  the  entire  pore  network  are  also  shown  for 
comparison.  In  order  to  convert  the  capillary  pressure  data 
obtained  with  octane  to  the  equivalent  water-air  capillary 
pressure  it  must  be  recalled  that  in  the  MSP  test  the  capillary 
pressures  are  measured  in  the  perfectly  wettable  standard, 
which  is  assumed  to  be  in  capillary  equilibrium  with  the 
GDL.  This  means  that  conversion  between  capillary  pres¬ 
sures  obtained  with  octane  and  water  is  done  by  means  of 
Eq.  (2)  with  both  contact  angles  taken  as  that  in  the  standard, 
which  is  zero. 

In  most  cases,  the  shapes  of  the  hydrophilic  pc-Snwp  curves 
are  quite  similar  to  their  respective  overall  /jc-vnWp  curves. 
This  supports  the  common  assumption  that  the  hydrophilic 
and  hydrophobic  pore  size  distributions  have  similar  shapes. 
In  cases  where  this  assumption  is  true,  a  single  overall  pore 
size  distribution  can  be  used  to  describe  both  pore  types  given 


that  the  relative  proportion  of  each  pore  type  is  known.  How¬ 
ever,  this  generalization  cannot  be  made  in  all  cases,  and  the 
Lyflex  felts  and  SGL  10BB  are  exceptions.  Although  the  dis¬ 
tributions  for  the  felt  materials  are  actually  similar  in  shape, 
some  differences  appear  in  certain  regions  of  the  curve.  This 
suggests  that  the  hydrophobicity  is  not  distributed  evenly  over 
the  entire  pore  size  range,  likely  due  to  the  presence  of  carbon 
powder.  The  plots  for  SGL  10BB  present  a  more  interesting 
result.  They  illustrate  the  similarities  of  the  10BA  and  the 
10BB  samples,  which  are  based  on  the  same  substrate  mate¬ 
rial.  The  overall  /;c-vnwp  curves  of  these  two  materials  diverge 
at  higher  capillary  pressures,  because  the  effect  of  the  micro- 
porous  layer  in  the  10BB  sample  is  superimposed  on  that  of 
the  substrate  material  which  is  common  to  both  GDLs.  It  can 
be  seen  in  Fig.  5  that  the  hydrophilic  portions  of  10BA  and 
10BB  are  virtually  identical  in  shape.  This  result  suggests  that 
capillary  pressure  curves  for  the  individual  components,  i.e., 
substrate  and  MPL,  can  be  combined  to  generate  the  overall 
capillary  pressure  curve.  This  possibility  is  investigated  in 
Section  3.2.6. 

3.2.4.  Correlation  for  hydrophilic  pore  distribution 

The  purpose  of  wet-proofing  the  GDL  is  to  render  a  por¬ 
tion  of  the  pores  hydrophobic  in  order  to  ensure  open  pores 
for  gas  flow  particularly  when  flooding  occurs.  This  treat¬ 
ment  confines  liquid  water  to  a  subset  of  the  overall  pore 
network,  which  is  of  great  interest  when  modeling  liquid 
water  behavior  in  GDLs.  Fitting  of  the  Brooks-Corey  and 
van  Genuchten  models  to  each  of  the  unimodal  hydrophilic 
distributions  yields  the  parameters  in  Table  5. 
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Table  5 


Model  fitting  parameters  for  unimodal  hydrophilic  pc-.Vwp  curves 


Material 

Brooks-Corey 

van  Genuchten 

Pcb  (bar) 

X 

RMS 

Pcb  (bar) 

m 

n 

RMS 

SGL10BA 

0.0636 

1.61 

0.38 

0.0918 

0.740 

3.843 

0.32 

SGL10BB 

0.387 

1.17 

0.33 

0.0580 

0.627 

2.683 

0.33 

Toray  090 

0.120 

4.06 

1.8 

0.1457 

0.8667 

7.499 

0.26 

E-Tek  Cloth  A 

0.0565 

0.950 

1.9 

0.0781 

0.5407 

2.177 

0.62 

3.2.5.  Leverett  J-function  correlation 

In  order  to  directly  compare  pc—sWp  data  from  different 
materials  it  is  necessary  to  account  for  the  effects  of  different 
material  and  fluid  properties.  Leverett  [32]  has  shown  that 
capillary  pressure  data  of  materials  with  similar  structure  can 
be  correlated  in  terms  of  the  /-function,  defined  as  follows: 

/(W=P"(^)1/2  (16) 

Since  most  GDLs  have  a  similar  carbon  fiber  structure,  pre¬ 
sumably  they  have  similar  pore  structures  so  that  the  Leverett 
correlation  should  be  able  to  provide  a  general  description  of 
the  Pc-s'wp  data.  It  is  not  necessary  to  include  the  contact  angle 
in  this  correlation  since  the  samples  are  all  made  of  similar 
material  and  the  capillary  pressures  have  already  been  con¬ 
verted  to  an  air-water  basis.  The  present  goal  is  to  determine 
whether  the  capillary  pressure  data  for  all  the  GDL  samples 
can  be  correlated  by  a  single  /-.vI1Wp  plot. 

To  convert  the  pc-sw p  data  to  the  /(.vnwp)  form,  the  total 
porosities  determined  by  MSP  with  octane  were  used  (see 
Table  3).  The  permeabilities  used  to  determine  /  are  listed  in 
Table  1 ,  while  the  surface  tension  values  for  octane  and  water 
were  taken  as  0.02175  and  0.0725  N  m-1,  respectively.  Any 
curves  that  showed  multimodal  distributions  were  excluded. 
The  plot  in  Fig.  6  shows  that  the  data  for  the  different  mate¬ 
rials  fall  reasonably  well  along  a  single  curve,  which  can  be 
fitted  by  VG-  or  BC-type  expressions.  The  fitting  parameters 
obtained  for  each  model  are  listed  in  Table  6.  The  parameter 


J  =  P(/(o(e/K)in) 


Fig.  6.  Comparison  of  experimental  /(s)-swp  curves  for  the  overall  and 
hydrophilic  pore  networks  with  Brooks-Corey  and  van  Genuchten  models 
for  all  unimodal  GDLs  tested. 


/c  was  introduced  into  the  BC  and  VG  models  in  place  of  pc b 
and  denotes  the  critical  curvature  of  the  /-function. 

The  success  of  the  /-function  correlation  is  particularly 
encouraging  given  the  diversity  of  the  materials  studied.  SGL 
10BA  and  Toray  090  have  similar  structures  composed  of 
randomly  oriented  carbon  fibers,  with  the  Toray  material  con¬ 
taining  no  PTFE  and  SGL  10BA  containing  5%  PTFE.  This 
indicates  that  the  Leverett  function  effectively  accounts  for 
the  differences  in  porosity  and  permeability  caused  by  the  use 
of  PTFE.  The  success  with  which  this  correlation  describes 
the  behavior  of  the  E-Tek  cloth  is  also  noteworthy.  The  woven 
substrate  is  obviously  very  different  from  the  randomly  ori¬ 
ented  SGL  and  Toray  papers,  yet  it  appears  that  this  difference 
does  not  affect  the  capillary  behavior  at  a  microscopic  level. 

3.2.6.  Microporous  layer 

It  is  interesting  that  the  contribution  of  the  MPL  to  the  pore 
size  distribution  of  the  SGL  1 0BB  sample  appears  to  be  distin¬ 
guishable  from  that  of  the  substrate.  This  can  be  clearly  seen 
by  a  comparison  of  the  pore  size  histograms  for  SGL  10BA 
and  10BB  compared  in  Fig.  7a  and  b.  In  the  range  between  10 
and  lOOnm,  the  10BB  sample  shows  additional  pore  volume 
due  to  the  microporous  layer.  In  Fig.  7c  the  capillary  pressure 
curves  obtained  by  MSP  with  water  for  both  10BA  and  10BB 
are  compared.  It  is  apparent  that  the  addition  of  the  hydropho¬ 
bic  microporous  layer  to  one  side  of  the  10BA  material  (i.e. 
to  create  10BB)  does  not  alter  the  /jc— .swp  behavior  of  the 
substrate  appreciably.  Based  on  this  information  it  should  be 
possible  to  determine  the  /)c-.vWp  curve  for  the  MPL.  How¬ 
ever,  this  is  not  a  trivial  matter  because  the  combined  effects 
of  the  added  MPL  material,  additional  pore  volume  in  the 
MPL  and  reduced  pore  volume  in  the  substrate  due  to  MPL 
penetration,  causes  SGL  10BB  to  have  a  larger  bulk  volume 
and  lower  overall  porosity  than  SGL  10BA.  Instead,  the  pore 
distribution  of  the  MPL  can  be  found  by  assuming  that  two 
unimodal  distributions,  one  for  the  substrate  and  one  for  the 
MPL,  can  be  added  together.  Thus,  the  saturation  .v®®  for  the 
10BB  GDL  can  be  obtained  from  the  following: 

swp(Pc)  —  +  (1  —  <p)swp  (17) 


Table  6 

Model  fitting  parameters  for  Leverett  J(s)  curve 


Brooks-Corey 

van  Genuchten 

Jc  A, 

RMS 

Jc  m  n 

RMS 

J(s)  0.6274  2.528 

0.44 

0.6982  0.7114  3.465 

0.51 
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Pore  Radius  (log  r)  [nm] 


Fig.  7.  Comparison  of  porosimetry  data  for  SGL  10BA  and  10BB.  (a)  Pore  volume  vs.  pore  size  histogram  of  SGL10BB  and  (b)  pore  volume  vs.  pore  size 
histogram  of  SGL  10BA.  Both  are  based  on  MSP  with  octane,  (c)  pc-sWp  curves  of  hydrophilic  pore  network  in  SGL  10BA  and  10BB  showing  similarity  of 
substrate. 


where  (p  is  the  proportion  of  the  total  pore  volume  in  the 
portion  of  SGL  10BB  not  intruded  by  the  MPL  material,  ,yBA 
the  saturation  of  SGL  10BA,  and  ,vBB  the  saturation  of  SGL 
1 OBB .  Values  of  tp  and  L  can  be  determined  by  minimizing 
the  difference  between  the  left  and  right  hand  sides  of  Eq. 
(17)  using  least  squares.  This  regression  yields  the  VG  model 
parameters  for  the  MPL  layer  as  well  as  (p  (Table  7).  As  shown 
in  Fig.  8,  this  approach  provides  an  excellent  fit  for  10BB  and 
Eq.  (17)  expresses  the  pc-^wp  dependence  for  this  GDL  very 
well. 

Estimates  for  the  penetration  depth  of  the  MPL  into  the 
substrate  and  the  porosity  of  the  MPL  can  be  found  by  inter¬ 
preting  the  <p  factor: 


q,  = 


vbb 

p 


Vs 

vp 

vB  eT 


(18) 


where  Vp  is  the  pore  volume  of  the  free  substrate,  Vp  B  the 
total  pore  volume,  Vg  B  the  bulk  volume  and  eBB  the  overall 
porosity,  all  pertaining  to  the  combined  GDL.  In  the  follow¬ 
ing  analysis,  it  is  assumed  that  the  portion  of  the  substrate 
penetrated  by  the  MPL  has  the  same  porosity  as  the  MPL. 
The  pore  volume  of  the  substrate  per  unit  bulk  volume  of 


Table  7 

Model  fitting  parameters  for  MPL  pc— ^wp  curve 


Material 

van  Genuchten 

Additional 

Pcb  (bar)  m 

n 

V 

RMS 

MPL 

8.443  0.5477 

2.211 

0.7938 

0.32 

substrate  can  be  taken  as  the  porosity  of  the  bare  paper  since 
the  porosity  distribution  of  the  substrate  is  not  altered  by  the 
addition  of  the  MPL  (see  Fig.  7a),  i.e. : 


Vs 

P  _  ,ba 

eT 


V, 


B 


(19) 


Eqs.  (18)  and  (19)  can  be  combined  to  eliminate  Vp,  yielding: 


Vj  <p4b 

t/BB  PBA 

VB  £T 


(20) 


Inserting  numerical  values  for  (p,  eBB  and  eBA  from  this  study 
into  Eq.  (20)  gives  a  substrate  bulk  volume  fraction  of  0.754. 


Air-Water  Capillary  Pressure  [Bar] 


Fig.  8.  Comparison  of  Eq.  (17)  with  experimental  Pc-^wp  data  for  SGL 
10BB.  The  curves  for  the  individual  MPL  and  substrate  contributions  are 
also  shown,  scaled  by  (p . 
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Fig.  9.  SEM  image  of  GDL  cross-section  obtained  by  freeze  fracturing  in  liquid  nitrogen,  (a)  SGL  10BB  showing  MPL  thickness  and  (b)  SGL  10BA  for 
comparison. 


The  bulk  volume,  Vg ,  of  the  substrate  per  unit  area  A  normal 
to  the  flow  direction  through  the  GDL  can  be  found  from 
knowledge  of  the  GDL  thickness,  i.e.: 

Vs  Vs  mpbb 

VB  _  VB  ,BB  _  .BB  ^ , 

A  ~  yBB'  -  PBA  1 

A  VB  ST 

where  the  SGL  10BB  thickness  fBB  is  reported  by  the  man¬ 
ufacturer  to  be  420  |Jim  (Table  1).  From  knowledge  of  the 
substrate  bulk  volume  fraction,  the  free  substrate  thick¬ 
ness  is  determined  to  be  317  p,m  and  the  MPL  thickness  to 
be  103  |im.  Assuming  the  quoted  SGL  10BA  thickness  of 
380  p,m  is  the  same  as  that  of  SGL  10BB  prior  to  the  MPL 
incorporation,  the  penetration  depth  of  the  MPL  into  the  sub¬ 
strate  is  estimated  to  be  63  |xm.  This  estimate  agrees  with 
visual  inspection  of  the  MPL  (see  Fig.  9).  Finally,  this  infor¬ 
mation  allows  a  value  of  0.72  to  be  determined  for  the  MPL 
porosity  f;^',pp  from  the  following  expression: 

fBB 

£MPL  =  (1_^eBB_  (22) 


This  value  for  ?;^1PP  confirms  the  expected  result  that  the  MPL 
has  lower  porosity  than  the  substrate. 


4.  Conclusions 

This  study  has  demonstrated  the  use  of  an  alternative 
method  of  porosimetry  capable  of  measuring  the  distribu¬ 
tion  of  the  total  and  hydrophilic  pore  volume  by  pore  size 
in  wet-proofed  GDL  materials.  The  total  porosity  distribu¬ 
tion  obtained  with  this  method  was  found  to  compare  very 
favorably  with  that  obtained  from  the  conventional  MIP  tech¬ 
nique.  The  experimental  /?c-.vwp  curves  of  the  overall  and 
hydrophilic  networks  for  GDL  samples  exhibiting  unimodal 
behavior  could  be  well  described  by  the  van  Genuchten 
model  over  the  entire  range  of  saturations.  The  shapes  of 
both  hydrophilic  and  overall  /;c-vWp  curves  of  these  GDLs 
were  also  similar  and  could  be  normalized  using  the  Lev¬ 
ered  function.  The  entire  data  set  for  these  samples  could  be 
correlated  using  the  van  Genuchten  model  to  provide  a  gen- 
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eral  equation  from  which  the  pc-sKp  curves  for  other  GDL 
materials  could  be  predicted.  However,  materials  with  car¬ 
bon  powder  or  a  highly  hydrophobic  layer  such  as  an  MPL, 
did  not  exhibit  this  behavior.  Separate  analysis  of  the  microp- 
orous  layer  of  one  GDL  sample  showed  that  its  behavior  could 
be  independently  described  by  the  van  Genuchten  model. 
Then,  this  could  be  superimposed  with  the  behavior  of  the 
standard  GDL  substrate  to  predict  the  /?c-.vWp  curve  of  the 
combined  substrate-MPL  GDL.  This  analysis  allowed  sev¬ 
eral  important  parameters  of  the  GDL  structure  such  as  the 
MPL  porosity  and  its  penetration  depth  into  the  substrate  to 
be  estimated.  The  van  Genuchten  model  that  was  used  to 
describe  the  pc-s Wp  curves  also  allows  the  relative  perme¬ 
abilities  to  be  estimated.  This  work  has  proposed  alternative 
relationships  for  the  modeling  of  liquid  transport  in  the  GDL 
of  PEMFCs.  It  is  anticipated  that  this  relationship  will  lead 
to  more  accurate  prediction  of  behavior  in  the  concentration 
overpotential  region. 
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